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We theoretically show that buckled two-dimensional graphene-like materials (silicene and ger- 
manene) subjected to a femtosecond strong optical pulse can be controlled by the optical field 
component normal to their plane. In such strong fields, these materials are predicted to exhibit 
non-reciprocal reflection, optical rectification and generation of electric currents both parallel and 
normal to the in-plane field direction. Reversibility of the conduction band population is also field- 
and carrier-envelope phase controllable. There is a net charge transfer along the material plane 
that is also dependent on the normal field component. Thus a graphene-like buckled material be¬ 
haves analogously to a field-effect transistor controlled and driven by the electric field of light with 
subcycle (femtosecond) speed. 


I. INTRODUCTION 

Novel Dirac materials such as silicene or germanene 
m are monolayers of silicon or germanium with hexag¬ 
onal lattice structures where charge carriers at the Fermi 
surface are, as in graphene, Dirac fermions EM- Re¬ 
cently, silicene has shown m promise for applications in 
electronics such as field-effect transistors (FETs) (2011231 
where, being a semiconductor, it has a natural advantage 
over graphene that is a semimetal. Below we will con¬ 
sider silicene but all qualitative results are also valid for 
germanene. 

In this paper we theoretically predict that a single 
monolayer of silicene (germanene) is controllable at op¬ 
tical frequencies by a normal component of the incident 
optical field just like the gate voltage controls channel 
current in FET. The main difference between silicene and 
graphene is that due to a larger radius of a Si (or, Ge 
in germanene) atom compared to a C atom, the corre¬ 
sponding hexagon lattice in silicene has buckled structure 
[ 23 j consisting of two sublattices that are displaced ver¬ 
tically by a finite distance L z ~ 0.5 A - see Fig. (TJa). As 
a result, silicene has large spin-orbit interaction, which 
opens up band gaps at the Dirac points (A so ss 1.55 — 7.9 
meV for silicene [H [55| and A so « 24 — 93 meV for ger¬ 
manene mm)- For graphene, the corresponding spin- 
orbit-induced gap is very small, 25 ^4eV [26). The buckled 
structure of silicene/germanene lattice allows also for the 
band gap to be controlled by an applied perpendicular 
electric field 1271: the band gap increases almost linearly 
with this electric field. 

Phenomena in silicene in a strong optical pulse field are 
illustrated in Figs. [l](b)-(e). A strong optical field causes 
electron transfer in the direction of the force :28, 29]. 
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In fact, a strong optical field in the ^-direction (normal 
to the silicene plane) decreases symmetry of the system 
from honeycomb (six-order, centrosymmetric) to trian¬ 
gular (third-order, non-centrosymmetric). This leads to 
appearance of effects such as optical rectification and in¬ 
duction of currents normal to the in-plane component of 
the applied electric field. 

Microscopically, the z component of the strong field 
causes transfer of electrons between the sublattices. As¬ 
sume for certainty that, for the chosen pulse, electrons 
are transferred from A to B. (Note that the change of the 
maximum field to the opposite, i.e., change of the carrier- 
envelope phase of the pulse by i r, would obviously cause 
an opposite transfer.) In the case of in-plane field F 2 D 
polarized in the y-direction, there is an electron transfer 
in both the y- and x-directions - see Fig. m- The sym¬ 
metry of the system dictates that with the reversal of F 213 
(for the same 2 -component, F z ) the y-current changes to 
the opposite but the x-current does not change, as shown 
in Fig.Qc). This implies, in particular, that the system 
causes optical rectification in the x-direction, which is 
due to the absence of symmetry with respect to the re¬ 
flection in yz-plane for either sublattices. 

Fundamentally different scenario takes place for F 2 d 
in the x direction - see Figs. m and (e). In this case, 
there is no current in the y-direction due to symmetry 
with respect to reflection in the xz-plane. With respect 
to field F 2 _d changing to the opposite, the x-current does 
not have any definite parity, which is rectification in the 
x-direction. 

To provide for the field-effect control of optical phe¬ 
nomena in silicene, the ^-component of the pulse electric 
field should be strong enough: F z > hu)/(eL z ) ~ 2 V/A, 
where oj is the optical frequency. Then, necessarily, the 
pulse should be very short, on the femtosecond scale, 
to allow the processes to be complete before significant 
damage to the lattice may have occurred - see Sec. |TT] 
below. For such fields, there may be partial adiabaticity 
(reversibility) set on, which we will show below in Sec. 
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FIG. 1: (a) Hexagonal lattice structure of 2D silicene. The 
lattice consists of two inequivalent sublattices labeled by “A” 
and “B”. Sublattices A and B are shifted in the 2 direction 
by distance L z . The angle of incidence of the pulse is 8. (b) 
Schematic of in-plane electron transfer induced by in-plane 
pulse electric field F 2 D directed along the y axis as shown. 
The curved red arrows indicate the electron transfer between 
the sublattices, (c) The same as (b) but for the opposite 
F 2 d- (d) The same as (b) but for the field directed along the 
x axis, (e) The same as (d) but for the opposite F 2 .d. The 
2 components of the pulse field has the same direction in all 
cases. 

II. MODEL AND MAIN EQUATIONS 

At the present time, record-setting ultrashort optical 
pulses have duration « 1.5 optical period .25] [3DJ, with 
duration of just a few femtoseconds. We will idealize 
and simulate such an ultrashort pulse with the following 
single-oscillation waveform, 

F(t) = F 0 e-“ 2 (1 - 2u 2 ) , (1) 

where F$ is the amplitude, which is related to the pulse 
power, V = cFq/Ak, c is speed of light, u = t/r , and 
r is the pulse length, which is set r = 1 fs. Note that 
this waveform has zero area, J^° F(t)dt = 0, which is 
required for a pulse propagating in far-held zone. 

We consider a p-polarized laser pulse with polarization 
direction parallel to the plane of incidence, orientation 
of which is determined by an angle ip measured relative 
to axis x. Here the xy coordinate system is introduced 
in the plane of silicene/germanene, oriented as shown in 


Fig. 0 a). The angle of incidence of the laser pulse is 
denoted as 8. 

Similar to graphene, the silicene/germanene mono- 
layer has honeycomb lattice structure, which is shown 
in Fig. Ha). The lattice has two sublattices, labeled 
“A” and “B”, and is determined by two lattice vectors, 
ai = a/ 2(-\/3,1) and a 2 = a/2(VS, — 1), where a is lat¬ 
tice constant, which is 3.866 A for silicene and 4.063 A for 
germanene. The distance between the nearest neighbors 
is a/\/3. The first Brillouin zone of the reciprocal lat¬ 
tice is a hexagon and is shown in Fig. ib). The points 
K = (2 7 r/a)(l/v / 3,l/3) and K' = (2 7 r/a)(l/v / 3,-1/3) 
are the Dirac points. In the buckled structure [see Fig. 
0a)], the 2 -shift distance is L z = 0.46 A and L z = 0.66 
A for silicene and germanene, respectively EUHg. 

For a graphene monolayer, where spin-orbit coupling 
is extremely small ( ss 0.03 meV), energy gaps at the 
Dirac points are correspondingly very small and can 
be set as zero for any practical purposes. Then the 
low energy spectra near the Dirac points are well de¬ 
scribed by the Dirac massless relativistic equation. For 
a silicene/germanene system, finite spin-orbit interaction 
opens up a much larger gap ~ 10 — 100 meV [27:. Such 
a gap in the energy spectrum of silicene modifies low- 
energy electron transport and interaction between elec¬ 
trons in weak magnetic fields 133j . However, this spin- 
orbit interaction is too weak and, in our case, can be 
safely neglected compared to characteristic energy scale, 
eF z L z , introduced by the strong electric field of the op¬ 
tical pulse in the buckled Dirac materials. At the same 
time, the buckled structure of a silicene monolayer in¬ 
troduces strong sensitivity of the system to the external 
normal field, F z m ■ Hence, based on this consideration, 
below in this article we disregard the spin-orbit interac¬ 
tion but take into account the buckled structure bringing 
about the sensitivity to the normal optical electric field. 

The Hamiltonian of an electron in silicene in the field 
of an optical pulse has the form 

H = H 0 + eF 2 d (t)v + (1 ^ y (2) 

where Ho is the field-free electron Hamiltonian, r = (a;, y) 
is a two dimensional vector, F 2d = (F x (t), F y (t)) = 
F(t) sin0(cos(/,sin</>), and F z (t) = F(t)cos9. Here the 
matrix form of the Hamiltonian corresponds to pseudo¬ 
spin, i.e., two components of the wave function ip a and 
ipB, which describe the amplitudes for an electron to be 
on the lattice site A and B, respectively. 

The field-free electron Hamiltonian, Ho, describes the 
nearest neighbor tight-binding model of silicene without 
spin-orbit terms. This Hamiltonian is exactly the same 
as the free-held Hamiltonian of graphene [57H5T1 and de¬ 
scribes the tight-binding coupling between two sublat¬ 
tices A and B see Fig. 0a). In the reciprocal space, 
the Hamiltonian Hq is a 2 x 2 matrix of the form [MJ [35] 

y - ( 0 7/(k) ^ rr , 

H °~\7 /*(k) ° )’ ^ 
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FIG. 2: (a) Hexagonal lattice structure of 2D sil- 

icene/germanene. The graphene lattice consists of two in¬ 
equivalent sublattices labeled by A and B. The vectors ai = 
a/ 2(v / 3,1) and a .2 = a/2(v / 3, —1) are the direct lattice vectors 
of silicene/germanene. The nearest neighbor coupling, which 
is characterized by the hopping integral 7 , is also shown, (b) 
The first Brillouin zone of silicene/germanene. Points K and 
K' are two degenerate Dirac points, corresponding to two 
valleys of low energy spectrum of silicene/germanene. The 
blue arrow shows in-plane (a;t/-plane) component of the time- 
dependent electric field of the pulse. The in-plane field, F 2 d, 
is characterized by azimuthal angle ip. 


where the hopping integral 7 is —1.6 eV for silicene and 
— 1.3 for germanene m, and 


The electric field of the optical pulse generates both 
interband and intraband electron dynamics. The inter¬ 
band dynamics introduces coupling of the states of the 
CB and VB and results in redistribution of electrons be¬ 
tween the two bands. For dielectrics, such dynamics re¬ 
sults in its metallization, which manifests itself as a finite 
charge transfer through dielectrics and finite CB popula¬ 
tion after the pulse ends [251 32 S3 ■ 

In the reciprocal space, the intraband dynamics is de¬ 
scribed by the acceleration theorem [46] . 

»!= eF <‘>- < 8 > 

This acceleration theorem is universal and does not de¬ 
pend on the dispersion law. Therefore the intraband elec¬ 
tron dynamics is the same for both the VB and CB. The 
time-dependent wave vector kj^q, t) of an electron with 
initial wave vector q can be found by solving Eq. ([ 8 ]) as 

e f* 

kr(q,i)=q+^y F(t!)dfi. (9) 


The corresponding electron wave functions are the well- 
known Houston functions [47) . 


$iq } (r ,*) = * 




k T (q,t) 


re 


/- 00 dt i-®c.[ k T(q.tx)] 


( 10 ) 


/(k) = exp + 2exp (-i^) cos (|») . (4) 

The energy spectrum of Hamiltonian T-Lq consists of the 
conduction band (CB) (w*, or anti-bonding band) and 
the valence band (VB) (w, or bonding band) with energy 
dispersion E c ( k) = ~ 7 |/(k)| (CB) and E v ( k) = 7 |/(kj| 
(VB). The corresponding wave functions are 

h?<r>-9( e -V) (5) 

and 

^ )(r) = 7f (e 7 ^)’ (6) 

where /(k) = |/(k)|e^ fc . 

The characteristic electron-electron scattering time 
T e - e in silicene/germanene is expected to be similar to 
the corresponding time is graphene, which is ~ 10 — 100 
fs [55M5] . The duration of the pulse in our problem 
(r p ~ 4 fs) is t p < r e _ e . Therefore, it is not unreason¬ 
able to assume that the electron dynamics in the external 
electric field of the optical pulse is coherent and can be 
described by time-dependent Schrddinger equation 


where a = v (VB) or a = c (CB). 

Using the Houston functions as a basis, we express the 
general solution of time-dependent Schrodinger equation 
0 in the following form 

^q( r ^)= 51 (11) 


Solution (111 is parametrized by initial electron wave vec¬ 
tor q. Due to the universal intraband electron dynamics 
in the reciprocal space, the equations, which describe co¬ 
herent electron dynamics in the pulse field, become de¬ 
coupled, greatly simplifying the problem. 

Expansion coefficients /3 Qq satisfy the following system 
of differential equations 


d/3 C g(t ) 
dt 


F 2 d(t)Q q (f) + eF z (t)L z (t. q) Q 
' ^ PnqWi 

d[d vc ^{t) _ 

dt 


F 2 d (t)Q* (t) + eF z (t)L z (t,q) 

l ----- 


Pcq(t), 


( 12 ) 


(13) 


where function L z (t, q), which is given by the following 
expression, 


(7) 

where Hamiltonian 7 i. of Eq. ([2]) has an explicit time de¬ 
pendence. 


L z {t. q) = L z e~* I-°° dtl{Ecl k Uq^i)]--EUk r (q, tl )]| ( 14 ) 


is specific to the buckled structure of silicene. It deter¬ 
mines the interband coupling induced by the perpendicu¬ 
lar component of the pulse electric field. Vector function 
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Qq(£) is proportional to the in-plane interband dipole 
matrix element, 

Q q (t) = D[k T (q,t)]e“«' / '-- d ‘ l{i5t:lkT(q ’ tl)I_E ’' [kT(q ’ tl)]} , 

(15) 

where D(k) = (D x (k) 1 D y (k)) is the dipole matrix ele¬ 
ment between the states of the CB and VB with the same 
wave vector k, namely, 


D(k) = (< c) 


er 



(16) 


of the velocity operator are 
07 


-ycc \jvv 

* t. ' r 


V3h 


sm 


ak x 


ak x 

7s 


sm I + 0k J cos ^ 


( 21 ) 


,, cc -» yuv a 7 ( ak x \ . aky 

v y=~ v y =h cos { 2V3 + ^) sm ^ 


( 22 ) 


Substituting Eqs. © and © into Eq. fl6| ), we obtain 
explicitly, 


D x { k) = 


and 


Dy( k) = 


ea 

l +cos (^) 

cos {Wt) 

— 2 cos 

' aky \ 

J. 

2V3 

1 +4 cos 


““(a 

) + cos 

' aky \ 

J. 


(17) 


sm 


(^) 


sm 


{ 3 ak 

V 2\/3 ) 




1 + 4cos 


cos 


( 3ak x A 
V 2\/3 J 


+ COS 


(4 1 ) 


(18) 

System of equations fl2[)-([T3| describes the interband 


electron dynamics and determines the mixing of CB and 
VB states in the electric field of the pulse. For undoped 
silicene, all VB states are initially occupied and all CB 
states are empty. Then the initial condition for system 
Eqs. ([I2|)-([l3) is (/3„ q ,/3 cq ) = (1, 0), and the mixing of 
the states of different bands is characterized by time- 
dependent component |/3 cq (f)| 2 . We also define the time- 
dependent total CB population by the following expres¬ 
sion, 


V™ = —i 


2 ay 

Vsn 


and 


cos 


cos 


ak x 

7 f 

ak x 

7 f 


+ (f) k cos 


aky 


, ,cv ■ 2a 7 • 
Vy = -i—r- sm 
y n 


ak. 


ak„ 


(23) 


~ s[n ( + < ^ k ) cos o • ( 24 ) 


The interband matrix elements of the velocity operator, 
V x v and V™, are related to the interband dipole matrix 
elements, V™ = iD x (k)[E c (k) — E v (k)\ /h and V™ = 
iD y (k) [E c (k) — E v (k)\ /H [|8J. Within the nearest- 
neighbor tight binding model, silicene has electron-hole 
symmetry, which results in the relation V y c = —V™. 

Let us denote current in the i direction induced by in¬ 
plane field F 2d in the j direction as J,- ; , where i,j = x, y. 
Similarly we denote charge transferred after the pulse 
ends through the system as Qij. This is determined by 
an expression 


Qij 



(25) 


V c (t)=^|^ cq (t)| 2 , (19) 

q 

where the sum is over the first Brillouin zone. The CB 
population, J\f c (t), characterizes the electron dynamics 
in silicene and determines whether the dynamics for the 
entire system is reversible or not. Namely, the dynamics 
is reversible if, after the pulse ends, the CB population, 
which is the residual CB population, is small compared 
to the maximum CB population throughout the pulse. 

Polarization of the system in a time-dependent electric 
field also generates electric current, which can be calcu¬ 
lated in terms of the velocity operator from the following 
expression 

j j(t) = JE E E ^ iq w v r a2 /w:*), (20) 

q ct±=v,c Oi2=v,c 

where j = x , y, and V“ ia2 are matrix elements of the 
velocity operator Vj = With the known wave 

functions ©-© of the CB and VB, the matrix elements 


The current can be expressed in terms of polarization 
P(£) of the electron system as J(f) = dP(t)/dt. Then 
the transferred charge is determined by the residual po¬ 
larization of the system as Qij = P l3 (t —> oo), where we 
introduced tensor indices for P similarly to those for J l;j 
and Qij. The transferred charge is nonzero only due to 
irreversibility of electron dynamics in the optical pulse 
field. For completely reversible dynamics, when the sys¬ 
tem returns to its initial state after the pulse, the trans¬ 
ferred charge would be exactly zero. 

III. RESULTS AND DISCUSSION 

A. Band Population Dynamics in Strong Pulse 
Field 

The principal distinction of silicene from graphene 
is that the sublattices, A and B, are separated “verti¬ 
cally” (i.e., in the ^-direction) by an appreciable distance, 
L z « 0.5 A - see Fig. 0A)- The strong field of the op¬ 
tical pulse causes non-perturbative nonlinear changes in 
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FIG. 3: (a) Pulse waveform as given by Eq. 0 for F 0 = 
2 V/A. (b) For excitation pulse polarized in the yz plane, 
CB population M c is shown as a function of time for pulse 
amplitudes Fq indicated. Incidence angle 0 = ±80°. (c) The 
same as (b) but for the pulse polarized in the xz plane with 
the direction of the maximum field, F, shown in the inset; 
incidence angle 9 = 80°. (d) The same as (c) but for 8 = 
-80°. 

the material. Such phenomena are sensitive to the max¬ 
imum field of the pulse, which is amplitude F 0 . For our 
choice of pulse Eq. 0, the maximum of the carrier oscil¬ 
lation occurs at the maximum of the pulse envelope, i.e., 
the carrier-envelope phase (CEP) is zero - see Fig. 0 a) 
illustrating the pulse waveform. 

The CB population, A f c , calculated in accord with Eq. 
for pulse polarized in the yz plane is displayed in 
Fig. 0b) as a function of time t for different field am¬ 
plitudes and incidence angle 9 = ±80°. Note that be¬ 
cause silicene is symmetric with respect to reflection in 
the xz plane, the results for both 80° and —80° are iden¬ 
tical. Two most prominent features of this dynamics are: 
(i) dependence on the pulse amplitude is very nonlinear, 
and (ii) the residual (after the pulse end) populations, 
M res) are close to the maximum populations during the 
pulse. The latter property is similar to that in graphene 
[29]. However, it is in a sharp contrast to that in silica, 
cf. Refs. [45] (theory) and {S5J [55) (experiment) where 
the residual CB populations are relatively small. This 
large residual CB population for silicene suggests lack of 
adiabaticity, which is likely due to a relatively small dis¬ 
tance of the transfer between the two sublattices in the xy 
plane, L xy = a/ (2\/3) ~ 0.7 A, in this case. Note that 
the adiabaticity parameter is <5 = Huj/ (eF y L xy ). Adi¬ 
abaticity requires S <C 1 while, in our case, even at the 
strongest fields, the adiabatic parameter is not too small, 

<S> 1. 

The response for the case of the pulse polarized in the 


xz plane is displayed in Figs. 0c)-(d). In the stark con¬ 
trast to the case of the yz polarization considered above 
in the previous paragraph, here there is a dramatic dif¬ 
ference between 9 = 80° and 9 = —80°. This is due 
to the violation in the reflection symmetry induced by 
the z component of the maximum field. For the case 
illustrated, this field promotes transfer of electrons pre¬ 
dominantly toward the B sublattice - cf. Fig. 0a). 

For 9 = 80° as shown in Fig. 0c), the x component 
of the maximum field, F x < 0, promotes transfer of elec¬ 
trons from left to right (in the direction x > 0) according 
to their negative charge - cf. Fig. 0d). The distance of 
transfer is the same as in the case of the yz-polarized 
field L xy = a/ (2x/3) ~ 0.7 A and adiabaticity is vio¬ 
lated since S = Huj/ ( eF y L xy ) > 1. Correspondingly, the 

residual CB populations are again close to their 

corresponding maxima during the pulse. 

Dramatically different behavior takes place for the re¬ 
ciprocal incidence, 9 = —80° where the CB population ki¬ 
netics is displayed in Fig.0d). For relatively weak fields, 
Fq = 0.5 —1.5 V/A, the kinetics is essentially irreversible, 
where the maximum CB population is attained at the end 
of the excitation pulse, similar to the case of Fig. 0c) 
considered above in the previous paragraph. In a sharp 
contrast, for stronger fields, F 0 = 2 — 3 V/A, there is 
partial reversibility: at the end of the pulse the CB pop¬ 
ulation is reduced by a factor of « 2 with respect to its 
maximum. This is related to improved adiabaticity, i.e., 
decreased adiabaticity parameter, <5 = fno/ ( eF y L xx ) < 1 
where L xx = a/\/?> ~ 1.4 A is the horizontal transfer dis¬ 
tance, see Fig. 0a). This distance is twice longer than 
for the case of Figs. 0c)-(d) corresponding to the polar¬ 
izations in Figs. 0bRc). 

Note that the adiabaticity in the case of Fig. 0d) is 
incomplete; for comparison, in the case of silica (quartz) 
a nearly perfect adiabaticity has been predicted and 
observed [351 S2]. This high degree of adiabaticity is 
most certainly related to a wide band gap, A g (see also 
Ref. [33]) and to a significantly larger lattice constant, 
a ~ 5 A, in quartz. Both these factors determine adi¬ 
abaticity, which is pronounced when fkv/Ag <C 1 and 
Hoj/(eF 0 a ) <C 1. Thus one should not expect near-perfect 
adiabaticity in graphene (cf. Ref. [29] ), silicene, and ger- 
manene where A g is negligible, and a is relatively small. 

B. Ultrafast Currents Induced by Strong Pulse 

Electrical current is due to displacement of charges 
caused by the applied pulse field. For free classical elec¬ 
trons, this current is proportional to their mean velocity, 
i.e., to the integral of the field, often referred to as vector 
potential, 

A (t) = -c f F 2d {t’)dt' . (26) 

J — OO 

In contrast to free electrons, as we have argued above in 
Sec. [TTT0 the strong field acting on the electrons in crys- 
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FIG. 4: Current dynamics in silicene subjected to strong field 
pulse. The broken line displays the shape of the excitation 
pulse. The numbers labeling the curves are the correspond¬ 
ing field amplitudes Fq (V/A). (a) The ^-component of the 
current density, J x , as a function of time t for the excitation 
pulse polarized in the xz plane with the direction of the max¬ 
imum in-plane field, Fq x . in the negative x direction [same as 
in Fig. [3jc)]. (b) Same as in panel (a) but for the opposite 
direction of Fq x [as denoted as XX, corresponding to Fig. 
id)]- (c) The //-component of the current density, J y , as a 
function of time t for the excitation pulse polarized in the yz 
plane, (d) Same as in panel (c) but for the x component of 
the current density. 


tal lattice of silicene causes effective symmetry reduction 
from honeycomb to triangular and, in particular, depen¬ 
dence of the electron dynamics on the sign of the maxi¬ 
mum field - cf. Figs. 3](c) and (d). The observed partial 
adiabaticity is also due to the presence of the periodic 
lattice and defined by its period in the field direction. 

The effective reduction of symmetry to triangular 
(where there is no inversion center) caused by the strong 
normal (z) field component causes the currents in the sil¬ 
icene lattice to be highly anisotropic and non-reciprocal 
as we show below in this Section. Let us denote Jxx 
an x-component of the current density induced by the 
field polarized in the xz plane with the maximum in the 
negative x direction as shown in Fig. §c). Similarly, 
we denote Jxx the x component of the current density 
caused by the field with the maximum in the positive x 
direction as in the case of Fig. [3](d). Note that gener¬ 
ally Jxx A — Jxx ( as would have been the case for free 
electrons) due to the low, triangular effective symmetry. 

Similarly, we introduce current density Jyy as the y 
component of the current density induced by the yz po¬ 
larized pulse. Note that in this case, the presence of the 
xz-symmetry plane dictates that Jyy = — Jyy. Interest¬ 
ingly enough, the in-plane field in the y direction causes 
also a current in the x direction [cf. Figs, [l] (b) and (c)], 


(a) 1.3fs (b) r=-0.7fs 



fe (l/A) 1 1 fe(l/A) 


FIG. 5: Electron momentum distribution |/3 c k| 2 in the first 
Brillouin zone in the CB for pulse with maximum field Fo = 
2.5 V/A with the Y-polarization of the in-plane field (the 
maximum in-plane field directed along the positive y axis), 
(a) Distribution at the moment of time t = —1.3 fs corre¬ 
sponding to the maximum negative current Jyy- (b) Distri¬ 
bution at the moment of time t = —0.7 fs corresponding to 
the maximum positive oscillation of current Jyy. 


whose density we will denote as Jyx ■ Note that due to 
the symmetry, this current is invariant with respect to 
inversion in the xz plane, i.e., Jyx = Jyx- 

In Fig. [4} we plot the temporal behavior of the current 
density for the four independent cases of the pulse polar¬ 
ization and current direction, XX, XX, YY, and YX, 
as indicated in the panels; the currents in all other cases 
are either related to these cases by symmetry, as pre¬ 
sented in the previous two paragraphs, or equal zero as, 
e.g., Jxy and -lyy- For the XX case shown in Fig. |4][a), 
in the relatively weak fields, F 0 < 1 V/A, the current 
density, Jxx , obviously, qualitatively follows the vector 
potential, A(t) reaching (negative) maximum at approx¬ 
imately quarter oscillation period and turning to zero at 
the maximum field (t = 0). Kinetics Jxx{t) is approxi¬ 
mately antisymmetric with respect to point t = 0, which 
shows that this process is nearly time-reversible. 

However, at higher fields, the behavior in Fig. |4][a) be¬ 
comes nontrivial. The first manifestation of this behavior 
appears at Fo = 1.5 V/A where instead of a pronounced 
minimum (maximum negative current) there is a plateau, 
which turns to a maximum for F 0 > 2 V/A. We attribute 
this behavior to electrons that are compelled by the field 
force to drift in the reciprocal space across the K-point. 
We will discuss this behavior in more detail in conjunc¬ 
tion with Fig. Qc) - see below. 

A phenomenon of fundamental importance is the loss 
of adiabaticity in higher fields, which manifests itself in 
the lack of anti-symmetry with respect to point t = 
0 in Fig. |4][a). Note that non-adiabaticity also im¬ 
plies irreversibility |53j and, consequently, violation of 
time-reversal symmetry (called also T-invariance or T- 
symmetry). This violation of adiabaticity is related to a 
gradual transfer of population between the A and B sub¬ 
lattices, as we discussed above in Sec. Ill A Such transfer 
is not instantaneous; one can estimate characteristic time 
it requires as t tr ~ TrH/(eL z F 0 ). For a high field used, 
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F 0 ~ 2 V/A, we obtain i tr ~ 1 fs. This is in a full quali¬ 
tative agreement with the results of Fig. |4^a) where the 
time-reversal asymmetry becomes pronounced for high 
fields and times longer than ~ 1 fs from the moment the 
pulse is applied. 

One of the consequences of the T-invariance violation 
are non-zero values of the transferred charge and of the 
residual polarization - see Eq. (25) and Fig. [6] and the 
corresponding discussion - violating the T-symmetry and 
adiabaticity. This implies that the system’s dynamics is 
irreversible (non-adiabatic), which may surprise one be¬ 
cause the system is completely Hamiltonian. This is due 
to the fact that the central frequency of the laser radi¬ 
ation, hui ss 1.5 eV, is close to the transition frequency 
between the electron states localized at the two sublat¬ 
tices, HAui ~ 7 rh/ttr = sL z Fq ~ 1.4 eV. This causes 
resonant absorption leading to dephasing - collisionless 
relaxation widely known as Landau damping Ell- 


Current kinetics for the XX case displayed in Fig.^b) 
is qualitatively similar to that for the XX case discussed 
above in the previous three paragraphs. However, the 
symmetry reduction caused by the nonlinear interaction 
with a controlled (zero in our case) CEP causes current 
■Jxx t° differ quantitatively from Jxx, which difference 
is pronounced in the second half-period (t > 0) where the 
T-asymmetry of the current becomes evident. The latter 
is due to the non-adiabaticity, already mentioned above 
in the discussion of Fig.|4](a): the transfer of the electrons 
between sublattices occurs during a finite period of time, 
ttr ~ 7r H/(eL z Fo) ~ 1 fs, comparable with half optical 
period in our case. 


The YY case illustrated in Fig. |4](c) is not related 
by crystal symmetry or other invariances to the XX 
and XX cases considered above. However, the kinet¬ 
ics of Jyy is qualitatively similar to, though quantita¬ 
tively different from, the previous two cases. Note that 
there is strict symmetry Jyy = —Jyy- Here also the 
T-symmetry is violated: the kinetics in the first and sec¬ 
ond half-periods is dramatically different. Note that in 
this case, current at the end of the pulse may not van¬ 
ish, which is certainly due to absence of collisions and 
other interactions in the model. Note that the electron- 
electron collisions is the fastest interaction-induced re¬ 
laxation process. However, it takes the electron-electron 
collisions ~ 10 — 20 fs in a similar two-dimensional sys¬ 
tem, graphene, to make an effect [41], which is too long 
time for our pulse whose entire duration is less than 4 fs. 

The results for current Jyx (in the x direction induced 
by the field in the y direction) are displayed in Fig. id). 
Note that exactly Jyx = Jyx due to symmetry. With¬ 
out an electric field applied, silicene is a center-symmetric 
solid. Therefore for low fields current Jyx should van¬ 
ish. This is, in fact, the case with a good accuracy for 
F 0 = 0.5 V/A, as one can see in the Fig. id). With 
held increasing, there is an increased current Jyx- Pre¬ 
dominantly, it is directed along the negative x axis, as is 
understandable from comparison with Figs. m and (c). 
Note that magnitude of this current is approximately an 



FIG. 6: Charge density Q transferred in the plane of the 
silicene layer as a function of the maximum pulse field, Fo- 
Four cases are displayed: Qxx, Qxx> Qyy , and Qyx, as 
labeled in the Figure. All other components of the charge 
density transferred are either zero, or related to these four 
cases. 


order of magnitude smaller than Jyy- 

The origin of the current oscillations for strong fields, 
F 0 > 2 V/A, in Figs.[S](a)-(c) can be understood from the 
electron momentum distribution. Consider for certainty 
the YY case, where the current is shown in Fig. |4|c). 
The corresponding momentum distribution for electrons 
in the CB, |/3 c k| 2 , for pulse held amplitude F 0 = 2.5 V/A 
is displayed in in Fig.[5](a) for moment of time t = — 1.3 fs, 
corresponding to the minimum (the maximum negative 
value) of the current, Jyy- At this instance, which just 
precedes the current oscillation, excess electron popu¬ 
lation (depicted by green) is concentrated at k y < 0, 
k x = 0. This excess population is formed due to held 
force eF y (t) > 0 that pushes the electrons across the K 
points into the second Brillouin zone in the extended zone 
picture; these electrons appear in the hrst Brillouin zone 
at the K point at k y « —1 A -1 , k x = 0. The second lo¬ 
calization of electrons is around a K 7 point at k x = 1 A -1 . 
This electron population is formed due to the lack of the 
center symmetry in the presence of strong held F z , i.e., 
it has the same origin at current Jyx described above in 
the previous paragraph. 

Dramatically different electron distribution is dis¬ 
played in Fig. §b) for t = —0.7 fs when current Jyy 
experiences the maximum upswing. This is caused by a 
signihcant number of electrons in the CB with k y > 0 
which appear due to a drift in the reciprocal space under 
force eF y > 0. These electrons make a positive contri¬ 
bution to the current (their group velocity v g < 0; cor¬ 
respondingly, due to e < 0, their contribution to Jyy is 
positive). 

The currents described above in conjunction with Fig. 
[4] cause transfer of charge across the system and accu¬ 
mulation of charges by the end of the pulse as given by 
Eq. ( p5| ). Such charge Q transferred through the system 
is displayed as a function of the field amplitude, Fo, in 
Fig.© for four independent combinations of the field and 
current directions, XX, XX, YY, and YX. A remark- 






able property of these results is that in all cases, except 
for YX, the charge transferred changes its sign as the 
field amplitude increases. This can be attributed to the 
increased number of electrons experiencing the Bragg re¬ 
flections at the Brillouin zone boundary, especially at the 
K-points, with the field increase. Thus this sign change of 
the transferred charge has the same origin as the current 
oscillations in Fig. [4] as described above. This charge ac¬ 
cumulated at the pulse end is an experimentally observ¬ 
able quantity just as the previous experiments on cur¬ 
rents in dielectrics [28l 152] . On the order of magnitude, 
this accumulated charge in Fig. [6] is Q ~ 1 fC//nn. For 
a ~ 1 pin focused spot, this gives a ~ 1 fC transferred 
charge. Such a charge is on the same order of magni¬ 
tude as in experiments Refs. |28J [55] and is, in principle, 
reliably observable. 

IV. CONCLUDING DISCUSSION 

The accumulation of charge Q transferred through the 
system implies a dramatic manifestation of symmetry vi¬ 
olation. This charge accumulation violates simultane¬ 
ously the parity symmetry (P-symmetry) and the charge- 
inversion symmetry (C-symmetry). This violation hap¬ 
pens due to the fact that our pulse is short and has a 
controlled CEP (zero in our case): the maximum field 
is reached at the maximum of the envelope (instance 
t = 0). Due to strong nonlinearity of the system for 
fields Fo > 1 V/A applied, this maximum field defines a 
selected direction in the system plane for the force act¬ 
ing on electrons. This causes the violation of the C- 
and P-symmetries. This is actually a general property 
for systems subjected to a short, strong, CEP-controlled 
pulses. It takes place in both two-dimensional solids such 
as graphene, silicene, germanene, and also conventional 
three-dimensional solids such as fused silica, sapphire, 
etc. In particular, it was a fundamental origin for the 
charge transfer in silica and quartz in original experi¬ 
ments [28] . 

A symmetry violation specific for silicene is related 
to the electron transfer between the sublattices caused 
by the normal field component F z , which effectively re¬ 
duces the system’s symmetry from hexagonal to trian¬ 
gular. This causes non-reciprocity: J\x 7^ Jxx and 
the appearance of a cross-current, ,Jxy A 0. Note that 
anisotropicity in the xy plane, Jxx A Jyy, is inherent 
in both silicene and graphene. 


Our zero-CEP pulse is T-symmetric; in the absence 
of the T-symmetry violation, the current should be T- 
odd, which would preclude the accumulation of charges 
after the pulse. However, we have seen from the results 
of Figs. |4](a)-(c) that the current is not anti-symmetric 
in time, i.e., there is a significant violation of the T- 
symmetry, which we attribute to the Landau damping. 
This is inherent in both graphene and silicene and is 
due to the absence of a significant band gap; this is 
in contrast to silica that is almost perfectly T-reversible 
l28l [30] . An additional contribution to T-irreversibility 
stems from the fact that the frequency associated with 
the electron transfer in the normal direction, eF z L z /H, is 
on the same order as the carrier frequency of the pulse. 
This causes resonant absorption of the excitation pulse 
and the Landau damping, specific for the silicene (and 
also germanene). If adiabaticity were present, it would 
have guaranteed reversibility and would have forbidden 
the charge accumulation. 

Finally, we note a close analogy of silicene with 
the field-effect transistor (FET) [50H22J. In FET, the 
gate field, applied normally to the conducting channel, 
changes the carrier populations in it and, thereby, con¬ 
trols its conductance. Analogously, in silicene, the nor¬ 
mal field component, F z . transfers carriers to one of the 
sublattices, A or B, thereby changing the system’s re¬ 
sponse to the in-plane field. A fundamental difference 
(and advantage) of silicene is that such a “device” works 
at optical frequencies, with the response time on the 
(sub)femtosecond scale. This opens a potential for many 
applications of silicene in future Petahertz-speed devices 
and applications. 
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